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Abstract: Interaction of charges in CCDs with the already accumulated charge distribution 
causes both a flux dependence of the point-spread function (an increase of observed size with flux, 
also known as the brighter/fatter effect) and pixel-to-pixel correlations of the Poissonian noise in 
flat fields. We describe these effects in the Dark Energy Camera (DECam) with charge dependent 
shifts of effective pixel borders, i.e. the Antilogus et al. (2014) model, which we fit to measure¬ 
ments of flat-field Poissonian noise correlations. The latter fall off approximately as a power-law 
^ with pixel separation r, are isotropic except for an asymmetry in the direct neighbors along 
rows and columns, are stable in time, and are weakly dependent on wavelength. They show varia¬ 
tions from chip to chip at the 20% level that correlate with the silicon resistivity. The charge shifts 
predicted by the model cause biased shape measurements, primarily due to their effect on bright 
stars, at levels exceeding weak lensing science requirements. We measure the flux dependence of 
star images and show that the effect can be mitigated by applying the reverse charge shifts at the 
pixel level during image processing. Differences in stellar size, however, remain significant due to 
residuals at larger distance from the centroid. 


Keywords: Photon detectors for UV, visible and IR photons (solid-state). 


Corresponding author. 






Contents 


m, Introduction 

^ Brighter/fatter effect in DECam 
g Model 

o A priori symmetries 

Assumed model symmetries 
[3.4| Fitting of coefficients 


3.2 Flat field covariances 


I 

I 

I 

I 

I 

0 

0 



1. Introduction 


An idealized felescope focal plane is covered by a grid of equal-sized pixels fhaf have a linear 
response fo fhe flux densify af fheir corresponding position in the sky. Charge counts contain 
Poissonian noise that is independent between any two pixels. The images produced by such a 
camera then, up to noise, represent the sky convolved with a kernel. The latter, called point-spread 
function (PSF), can be measured from the observed profiles of point sources. The only difference 
between point sources, e.g. stars, of different flux levels is an amplitude by which the PSF profile 
is re-scaled. 

Real CCDs with large well depth exhibit deviations from this picture, two of which have been 
known for some time. In exposures of homogeneous flux density (flat fields), the variance V at 
moderate to high flux levels is not linear in the flat level /i, unlike what is expected for a Poissonian 
process O- Secondly, the size of star images has been found to increase with flux (known as the 
brighter/fatter effect), e.g. in science verification data from the Dark Energy Camera (DECam), 
but also in many other cameras (e.g. MegaCam and ESST candidate sensors [Q], Euclid candidate 
sensors [26| and the Wendelstein Wide Eield Imager p^]). 

The brighter/fatter effect is particularly problematic for ongoing and future large weak leasing 


programs such as the Dark Energy Survey (DES, [plU), the Hyper Suprime-Cam Survey (HSC, 


[]24|]), the Kilo Degree Survey (KiDS, [^), or the Earge Synoptic Survey Telescope (ESST, Ql9|]). 
The reason for this is that a systematic misrepresentation of the PSE in bright stars (that are used 
for PSE modeling) biases galaxy shape measurement, to a level that can exceed the acceptable 
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systematics (cf. Section 3 and, e.g., |[T]]). Careful correction of this effect is therefore of great 
importance. 

Two physical phenomena have been proposed and examined with analytical and numerical 
models that could be the cause for both observations: (1) lateral electric fields due to 

differences in electrical potential between pixels of different charge count (put more simply, mutual 
deflection of accumulating charges) and (2) an increase in lateral diffusion due to lower drift fields 
in pixels fhaf are partially filled [|T^. 

iQ 1^ (hereafter A14) have proposed a phenomenological model for charge self-inleraclion 
in CCDs fhaf describes bofh fhe variance non-linearify and fhe brighler/falfer effecl. In fheir model, 
accumulating charge shifls fhe effective pixel borders. In a flal field image, a pixel fhaf has af 
any poinf during an exposure collecfed more charges fhan ifs neighbors due fo positive confribu- 
fions from Poissonian noise will, effectively, shrink. This decreases fhe variance and infroduces 
a positive covariance with neighboring pixels. The charges already present in the image of a star 
shift pixel borders inward, such that the measured image on the pixel grid increases in size. Using 
laboratory and on-sky measurements with several different cameras, A14 have shown this model 
to reproduce the observed effects reasonably well. Both lateral electric fields and increased laferal 
diffusion due fo lowered longifudinal elecfric fields are described by fhe same phenomenological 
model. 

One of fhe insfrumenfs sfudied by A14 was DECam |^, ^ |^, 21, 13, 10|, fhe 3 sq. deg. 

camera mounfed af fhe prime focus of fhe Blanco 4m felescope af fhe Cerro Tololo Infer-American 
Observafory (CTIO), currenfly used for DES and communify observing programs. The mosaic 
consisfs of 62 fully depleted 2Kx4K science CCDs wifh a fhickness of 250/im. During fhe science 
verification phase of DECam, we defecfed a flux dependence of observed sizes of sfar images, in 
line wifh fhe A14 sfudy, af a level relevanf for DES science requiremenfs. 

In fhis work, we aim af correcfing fhe effecls of charge self-inleraclion in DECam images af 
a level tolerable for DES weak tensing science. In Section 2, we qualifafively characterize fhe 
brighler/faller effecl from DECam sfar images. Secfion 3 infroduces fhe A14 model and explains 
our measuremenl of fhe model paramelers from flal fields. In Secfion 4 we discuss fhe effecl 
of charge self-inleraclion in DECam on galaxy shape measuremenls for weak tensing. Secfion 5 
describes fesls of our correcfion for charge self-inleraclion by reverse application of fhe model. We 
summarize our findings in Section 6. 


2. Brighter/fatter effect in DECam 


In fhis secfion, we make a phenomenological description of fhe flux dependence of fhe PSE in 
DECam dala. 

We illuslrafe the effect in 2D by comparing the profiles of brighl and fainl slars in DECam 


images. The lefl panel of Eig. 10 shows fhe residuals of normalized profiles of slars af a peak flux 
(SExtractor 0 parameter FLUX_MAX) of 5000 and 20000, i.e. (PSE 20000 — PSE 5 ooo)/PSE 5 ooo 
where PSE^? is fhe PSE model made from slars around FLUX_MAX = N. Elux is missing in fhe 
cenfral region of the profile and re-appears on an annulus af ss 2 — 4 pix radius. 

We use fhe flux deficif in fhe cenfral pixel of fhe normalized PSE model as a mefric fo answer 
a number of queslions. To fhis end, we use i band images of the globular cluster ft) Centauri taken 
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Figure 1. Flux deficit of normalized PSF model in central pixel relative to a PSF model made with stars 
in peak flux (FLUX_MAX) bin centered on 5000 ADU. We find that the flux deficit is approximately linear 
in flux level (solid lines show best linear fit, dashed lines 68 per cent confidence intervals), independent of 
exposure time (green squares: 50s, blue triangles: 80s), and increases when signal chain non-linearity is 
corrected for (dark green/blue: corTected, light green/blue: uncorrected). 


with DECam during science verification on Feb 02 2013 with homogeneous seeing conditions. 
We process these as described in Section |5|, without the reverse charge shift model applied. Fig. |] 
shows results as a function of FLUX_MAX for star images with two different exposure times and for 
a data reduction scheme applying or not applying, respectively, a correction for the known signal 
chain non-linearity of DECam. ^ We find that 

• the relative flux deficit in the center of bright stars due to the brighter/fatter effect is well 
described as linear in flux, as indicated by the solid lines that are good fits to the data points; 

• the amplitude of the flux deficit does not depend on exposure time, but only on the measured 
count level, as visible from the agreement between the blue and green lines (dashed lines 
give confidence limits for the linear fit); this indicates that the charges must be moved while 
being collected, not while they are waiting for readout in the pixel potential well; 

• the effect is amplified when the high level signal chain non-linearity, a property of the readout 
electronics, is corrected for; relative photometry of stars in frames with different exposure 
times shows that correction for high level signal chain non-linearity is necessary for restoring 
the physical charge level present in the CCDs (cf. Bernstein, presentation at PACCD2013 
conference^); we will therefore apply the signal chain non-linearity correction in all follow¬ 
ing analyses. 

*The non-linearity correction was determined by comparing the mean count levels of dome flats of different exposure 
time. At the 20k ADU level charge counts are approximately 1% lower than expected from linear response. 

^See https://indico.bnl.gov/event/cosmo2013. 
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Figure 2. Sketch of a DECam CCD and definition of our coordinate system, with the x and y axis parallel 
to the channel stops and the serial register, respectively. The split of the chip into two amplifier regions is 
indicated by the central horizontal line. 


3. Model 


We briefly introduce the model of A14 before we discuss our approach of constraining the model 
parameters in Sections ^ and 

The basic idea is that charges q^i in any pixel {k,l) influence the path of newly generated 
charges on their way into the pixel well.^ As a result, they induce a shift in the effective border of 
a pixel (/,y), i.e., they change the mapping of physical area on the chip onto the pixels. 

Consider pixel {ij) and its four borders, which we label as L = (—1,0), /? = (1,0) (left/right 
border to neighbors along x axis) and T = (0,1), B = (0,-1) (top/bottom border to neighbors 
along y axis). The shift in any of these we call 5f-, where X G {L,/?, T,B}. For a definition of our 
coordinate system, refer to Fig. ^ 

We assume the border shift to be linear in the surrounding charge (including charge in pixel 
(/, j) itself) and dependent on the lag {k — i,l — j) between pixel (t, j) and {k,l), 


kl 


(3.1) 


The parameters of the model are the shift coefficients afj. In all cases, 5^j shall be defined such that 
the positive direction is outward from pixel {i,j). 

The change in charges observed in pixel {i,j) is the sum of and the corresponding 

terms for L,T,B. Here, we have approximated the charge density at the pixel border by the mean 
of the neighboring pixels and defined the d as fractions of a full pixel width. The observed charge 
Qij then becomes 


QiJ ~ ‘lij + 


<lij + R , ^ij + ‘li-t.j „L 
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, + „ , ‘? 0 '+ 

+ -5- 2_i‘^k-i,l-ffikl + -X- ^<^k-i,l-ffikl 


(3.2) 


^Note that this could be due to a combination of lateral fields and changes in lateral diffusion. 


















3.1 A priori symmetries 

Assuming conservation of effective area, the area lost in one pixel is gained by its neighbor, 


xR — xi 

Oij - -Oi+l,j 

xT _ sB 
Oij - -Oij+\ 


(3.3) 


Eqn. 3.1 fulfills these equalities for general charge distributions q if and only if 


L, R 

Qij j , 

fjB — _„T 

‘^ij ■ 


(3.4) 


The number of free parameters is further reduced by two parity symmetries. The first is in¬ 
variance of charge shifts under change of sign of one coordinate axis. 


o,±i o.±i 

^ij =^-iJ ■■ 
±1,0 ±1,0 


(3.5) 


Another parity symmetry is that charges on opposite sides of the pixel border cause opposite shifts. 


o,±i _ _ o,±i 

±1,0 ±1,0 


(3.6) 


With these symmetries, the free parameters of the model are reduced to, e.g., afj for i >0, j>0 
and afj for i >0,j> 0. Note that we discuss further empirically assumed symmetries, necessary 


to unambiguously constrain the parameters, in Section 

3.2 Flat field covariances 


The mixing of independent Poissonian processes q into observed counts Q in Eqn. causes a 
correlation of noise in flat-field images. A14 collecf fhe ferms af firsf order in a fo find^ 


Cov(eoo,Go')=2At' £ 


(3.7) 


X=T.B,L,R 


assuming /r is fhe mean charge counf level and Poisson variance [Var{q) = /i] of charges in pixels in 
fhe flaf image. We can fhus use measuremenfs of flaf field covariances fo consfrain fhe a paramefers 
under suifable symmefry assumpfions. Nofe fhaf, in addifion, fhere is a change in variance, which 
can be wriffen using fhe above symmefries as 


AVar(eoo) =-4/1 (aio + ^oi), 


(3.8) 


equal wifh negafive sign fo fhe sum of all covariances infroduced befween Qqq and fhe surrounding 
pixels. 

Nofe fhaf if we converf Q and /r to quantities in ADU with some gain g, Eqns. and still 
hold with the same coefficients a. In Eqn. ^ we have to multiply the coefficients a ^ a' = ga to 
get the change in observed ADU. 


‘*Note that in their definition the a are only half as large, origin of the factor of 2 here. 
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Figure 3. Checks for systematics in measured covariances, shown for DECam chip N1. Left: median versus 
mean value of ~ 1200 measurements for each lag (red: amplifier A, green: amplifier B). Center: covariances 
measured at lag {i,j) versus Right: covariances measured in amplifier A region versus amplifier B 

region. From these tests we conclude that the effect of outliers in measured covariances, artifacts, and 
differences within the chip are negligible. 


Measurement pipeline 


We measure the covariances as follows, using all 10 s (« 15000 ADU) dome flat fields in the r 
band from the full first year of DECam observations (2013-08-15 until 2014-02-09). We perform 
overscan and bias subtraction based on the nightly median bias. Inter-CCD cross-talk between the 
two amplifiers and the signal chain non-linearity are corrected for. Since cosmic rays are systematic 
contaminants of the observed covariances, we run SEXTRACTOR on each flat image to detect all 
groups of three or more pixels 3 a above the background. Eor each night, we mask the union of 
these 2x2 box-convolved object masks, a bad pixel mask, and edge distortions [18, 22, ^]. We 


re-scale each flat image so as to make median count levels of all frames taken in the same night 
match, which changes the levels below the per-cent level. Erom each pixel of the flat images, we 
subtract the mean of that pixel over all n flats of the night. 

Erom these residual images, we estimate covariances out to lags of (8,8). We re-scale the 
residuals by ., / to correct the bias of the maximum likelihood ensemble covariance estimator. 

•' y n—l 

Since each chip has two readout channels, we make measurements separately for the 118 halves of 
chips in DECam,^ i.e. on a chip by chip and amplifier by amplifier basis. 

Covariances are predicted by the A14 model to scale with the square of the count level jj. 


(cf. Eqn. |3.7| ). We have verified this behavior of measured covariances from a PTC series of pairs 
of flat fields of different exposure times. Eor all following results, we therefore normalize our 
covariance measurements by to make them comparable despite somewhat varying flat levels 
between different nights and filters. 


Checks for systematics 

Each of the Ril200 flat images yields an independent estimate, such that we can perform a number 
of systematic checks. Eor each chip, we make comparisons between different measurements that 

^Of the 62 chips, two (N30 and S30, i.e. CCDNUM 61 and 02) are damaged and one (S7, CCDNUM 31) has a time- 
varying offset, which is why we exclude them from our analysis. 
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Figure 4. Scaling of covariance and measured gain with count level, as measured in a series of pairs of flat 
field exposures with different exposure time. Left panel; covariance at (1,0) lag, with f' scaling around 
the f = 10 s point indicated by red line. Right panel; measured gain g — fi/V, with the red line as a linear 
best fit to the data. 


should give the same result; 

• Since the statistics of our covariance estimator is very close to Gaussian, we expect no sig¬ 
nificant differences between using the mean and median of the « 1200 independent mea¬ 
surements. Outliers, however, would influence the mean more strongly than the median. 

• Chirally-symmetric lags (e.g. (1,1) and (1,-1)) measure covariances on disjunct sets of 
pairs of pixels that are, however, predicted to be the same from the model. 

• If chips are homogeneous and the model is correct, then Cov{Q oq, Qij )/should agree 
between the two halves of the chip read out by different amplifiers (cf. Fig. ^ despite their 
different gain values. 

We find all these measurements to be consistent within the errors with systematic effects that are 
small compared to our targeted accuracy, out to large lags (cf. Fig. 


Flux scaling 

Covariances of the flat field Poisson noise are predicted to scale with according to the model 
of Eqn. 3.7. Measured gain, i.e. the ratio of count level and variance, is predicted to rise linearly 
with flux due to the lowered variance. Other effects, such as correlated read noise, charge transfer 
inefficiency or other non-linearities in the signal chain, could cause different scalings. We test this 
by measuring covariances and gain at a range of count levels. 

To this end, we analyze a photon transfer curve (PTC) series of pairs of flat field images with 
exposure times of 0.05 ... 30s at a light level of ss 1500ADU s^^ Figure ^ shows covariances and 
gains measured from each pair, averaged over all valid chips. At very low levels, there is a floor 
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Figure 5. Left: Test for time stability of flat field covariances for DECam chip N1. Shown are measurements 
in flat frames of each month of the DES season of the covariance at lag (1,1) (red and green for amplifier 
A and B, respectively) and (1,-1) (magenta and blue). The four measurements are offset on the time axis 
for readability. No significant time variability is observed. Right: Wavelength dependence of pixel-to- 
pixel covariances as measured in DECam g, r, i, z, and Y flat field images. Measurements are plotted at 
the respective central wavelength for lags (1,0) (solid symbols), (1,1) (crosses) and (0,1) (open symbols). 
Intervals of ±10% around the r band measurements we use to fix the parameters of the model are indicated 
by dashed lines. 


of covariance below Cov < 0.01, which could be due to a small correlation of noise in the signal 
chain or a low-level non-linearity. At larger levels, correlations rise more quickly than predicted, 
potentially due to saturation and charge transfer effects. In the range of 2... 15 s exposure time, 
however, the Cov oc scaling describes the measurements well. 

Likewise, measured gain is predicted by the model as linearly increasing with count level. 
Gain measurements in flats of ^20 s exposure time overestimate the true gain by sslO % and we 
correct the gains used in the DES data management system ||^ for this effect. 

Wavelength dependence 

A14 found the noise correlations to be achromatic for the e2v-250 CCD run with nominal voltage 
configuration. This is in line with an effect that acts on the charges primarily in the last few /rm of 
their drift path and therefore is almost independent of conversion depth. We measure covariances 
on g, i, z, and Y band flats from 20 nights to test for a potential wavelength dependence. Results 
for the three innermost lags, averaged over all 59 chips, are shown in the right panel of Fig. ^ Both 
g and i measurements are consistent with the r band baseline model within 10%. For the (0,1) lag, 
the increase from g to Y band is at the 40% level between the g and Y band. In (1,0) and (1,1) 
there are indications of an increasing trend, although at a weaker level. 

Time stability 

We also test our measurements for time stability (cf. Fig left panel) by binning the measurements 












Figure 6. Pixel-to-pixel covariance as function of distance r = a/ f- + p. The long axis and readout direction 
is defined as x. Correlations with pixels along the same x or y coordinate are plotted in blue and red, 
respectively. A power-law with arbitrary amplitude is indicated by the dashed green line. 




■—r 

1 1 T 

^ ^ r 

n 




1 1 T 

^ ^ r 





1 1 T 

^ ^ r 





X 




1.2 



* 



1.8 




- 


7 








- 








- 










X 








- 

\ 










* 




i 




o 



K 

X 


o 




X X 

X X 

o 


* 




X 

6 





X 



5 “ 

* « 

I 

X 

1.6 




- 

o 















s 


. 



















•> 


X 

* 



> 






> 




X X 

■ 

(1 






o 






o 




1 

X 

o 





X 

o 






u 




1 

8 “ ■ 


5 

- K 


* « 



0.8 


X 







X 















1.4 









" 1 




■ i 







» * 

i 




1 

1 1 1 

1 1 1 

1 1 



• 

1 1 1 

1 1 1 

1 1 • 



• 

1 1 1 

1 1 1 

1 1 • 



1 

2 3 4 

5 6 7 

8 9 



1 

2 3 4 

5 6 7 

8 9 



1 

2 3 4 

5 6 7 

8 9 




lot 

number 





lot 

number 





lot 

number 



Figure 7. Cov/jX^ measured from flat fields for each of the DECam chips, plotted as a function to CCD 
production lot. Marked in red are the two high-resistivity lots. Left; (1,0) lag (left/right neighbor along 
readout direction), center; (0,1) lag (top/bottom neighbor across the channel stop), right; (1,1) lag (diagonal 
neighbor). 


on a monthly basis. We find no evidence for time variation in any of the 59 chips used in our study. 

Results 

We show covariances as a function of lag distance r = \JP- + p, averaged over the 59 chips, in 
Fig-i We measure positive correlations with a signal-to-noise ratio of ss 15 even at the outer¬ 
most point r = V128. Off-axis covariances show a smooth power-law drop-off with radius, with 
indications for a break at around r k, 5. The behavior of on-axis coefficients is slightly different, 
with a higher amplitude and a low outlier at the (0,1) lag, i.e. for neighbors across the channel 
stop. We note that the ratio Cov(l,0)/Cov(0,1) « 5 is in the range the [ |2^ simulations predict for 
reasonable channel stop barrier dipole moments. 
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We test for variations of the level of covariances between different chips. Since all chips 
in DECam use the same nominal bias and clock voltages, the reason for such differences would 
have to be connected to chip properties. The science chips used in DECam were produced in nine 
different production lots. Eor simplicity of notation, we label the two high resistivity (10.6-14.2 kfl 
cm) lots 107419 and 112094 as 1 and 2 and the remaining lots (with lower resistivity of 4.3-6.2 kfl 
cm) 123194, 123195, 124750, 124753, 135959, 135960 and 135961 as 3-9 (in this order). Eig. ^ 
shows results for the three nearest lags. We find variations in amplitude at the ps 20% level from 
lot to lot. Variations within the lots are significantly smaller, although with outliers, particularly at 
the (1,0) and (1,1) lag in lots 3 and 5. 

The two high resistivity lots have covariances « 20% higher (lower) than the low resistivity 
ones in the (1,1) (the 0,1) lag. Earger lags show a similar resistivity dependence as (1,1). Eor 
the (1,0) lag, no strong resistivity dependence is observed. This is evidence for a mixing of two 
effects with different dependence on resistivity and distance. While this suggests that resistivity is 
a relevant factor, other unknown lot-to-lot variations that coincide with resistivity by chance cannot 
be excluded, particularly since the high resistivity lots were cut from a different boule than the 
remaining CCDs. 


3.3 Assumed model symmetries 


The parameters of the model are the shift coefficients a from Eqn. |3.1| . They are connected to the 
measured covariances (see Section 3.2) by Eqn. With only the symmetries of Eqns. 3.4ff, how¬ 
ever, the shift coefficients are not unambiguously constrained by flat field covariances (cf. A14). 

This is most comprehensibly exemplified by the following degeneracy. Einding all covariance 
terms that contain any afj with i < 2 and j <2 


ACov(2oO)2oo) — 0+ '^0,1) 

Cov(2oo,Go,±i) = 2iJ.^{all — 2af i) + ... 

Cov(2oo,G±i,o) = 2At^(afo-2a[i)-h... 

Cov(2oo,2i,i) = 2/r^(af 1 -ha[i)-h... (3.9) 


we see that all covariances (including the variance defect ACov) are degenerate under the transfor¬ 
mation 


a^ 0 A 
afu-^ 
ah - A/2 

a[i-hA/2 (3.10) 

with an arbitrary constant A, which compensate the smaller shifts of one pixel border by larger 
shifts of the perpendicular ones. Clearly, however, the effect on inhomogeneous surface brightness 
images is not invariant under the above change on model coefficients. 

Our approach is to assume two additional symmetries (the rotational and projection symmetry 
described below) that are suggested by the covariance measurements. Eor the above degener¬ 
acy, this entails setting ah = af 1 and modelling the observed asymmetry in Cov(2oo! 2o.±i) / 
Cov(2oo,2±i,o) by differences in ah and ah only. 


T 

'^ 0,1 

ah —)• 
ah —)• 
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Rotational symmetry 

All off-axis lags rotated by 7 r /2 have covariances measured to be equal at the level of agreement 
between the A and B sides of chips, i.e. Cov{Qoo,Qij) = Cov{Qoo,Qji) for i > 0, j > 0. This 
symmetry is predicted by the model if 

afj = a]i'^i>0J>0. (3.11) 

If we assume this symmetry, we forfeit almost half of our observables (because the covariances at 
lag {i,j) and {j,i) are now two measurements on the same combination of shift coefficients). In 
order to break the degeneracies, it would suffice fo sef afi = ' > 0 , wifh a sef of assumed 

rafios c,. Due fo fhe agreemenf of measured covariances af rofafed lags, however, fhis would jusf 
yield a sfafisfically consisfenf (yef noisier) model. 


Projection symmetry 

The lateral electric field due to two charges at the same distance from a border should result from 
an electrostatic force that has the same amplitude but a different direction.^ The most similar 
pairs of charges in terms of distance are ones at lag {ij) and If we assume that the border 
displacement is proportional to the component of the force projected onto the normal vector of a 
border, we have (for / > 1 or 7 > 1 ) 


«7,i = , 

(3.12) 

where 


7-0.5 /(/-0.5)2+ 72 

/-0.5y (7-0.5)2 + /2 • 

(3.13) 

Eor a^, we write the analogous 


7 

II 

(3.14) 


As we find from a first iteration of fitting, the fall-off of shift coefficients with radius is close 
to a power-law with slope a = —2. We use that to refine Eqn. 3.13 as 


r{ij) = r{i,j) 


( /(7-0.5)2 + /2\“ 


7 -0.5 / (/-0.5)2 + 7-2 
/-0.5 V(7-0.5)2 + /2 J 


(3.15) 


and insert / instead of r in Eqns. p.l2| and p.l4| for our final model. 


3.4 Fitting of coefficients 

Assuming these symmetries, we have reduced the vector a of independent coefficients to af^, aj j, 
af J and afj for 1 < i > 7 > 0 (a total number of + 1) ~ \{n+ \){n + l) when going out 

to a maximum of n pixels distance). This matches the dimension of the vector c of independent 
measured covariances at the same lags. 

®Note that the situation is different for lateral diffusion, which is affected by the change in the longitudinal electric 
field. For distant enough pairs, however, lateral electric fields appear to be the dominant effect. 
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Out to lags of A = 3 we use the direct measurement of covariances of each chip individually. 
At larger distances, the chip-wise signal-to-noise ratio of the covariances is low. We therefore 
assume a power-law fall-off of covariances with radius, as observed in Section 

Cov(/,;)=A(/2 + /)-/5/2 (3.16) 


with j3 = 2.5. We fit three independent amplitudes (one for lags along the .r and y axis and one for 
off-axis pixels) to the covariances measured between 3 and 5 pixels separation. Covariances out to 
a maximum lag of 25 are extrapolated from this model. 

Using Eqn. ^ and the symmetries we express c as a linear function of a, c = Ma. We invert 
to find a = as a function of the measured and extrapolated vector of covariances. 

We note that this approach is somewhat different from the one proposed by [[I^, who model 
the radial fall-off of the shift coefficients in LSST candidate sensors by an exponential integral. 


Tests 


We perform two consistency tests of the fitted model. 

By comparing to measurements (Fig. left panel), we confirm that the model correctly repro¬ 
duces (at lags < 3) and extrapolates (at larger lags) the observed covariances. 

We further compare the fitted coefficients to our gain measurements from the PTC series 
(cf. Fig. The deviation of g from the true gain g is predicted to be linear in flux with a slope of 
4g(af Q -|-aQ j). We show the measured slope from the PTC series plotted against model parameters 
for all chips in Fig. ^ (right panel), finding the two to be in agreement at the 10% level. Note that 
the fitted + are quite sensitive to changes in the shift coefficients at large radii, i.e. the 


maximum lag used or the power-law coefficient of Eqn. 3.16 


Results 

We show the fitted shift coefficients of one exemplary chip in Fig. ^ The fact that for rows 
at Ay > 2 there is a maximum at Ar > 1 (or, equivalently, that curves for constant Ax cross as a 
function of Ay) is due to projection effects: at larger Ax the separation increases, but the electrostatic 
force is aligned closer to the normal of the pixel border. Using the projection angle 6, we correct 
for this effect by normalizing with cos G = [Ax — 0.5]/[y^(At— 0.5)^ + Ay'^]. 

The de-projected coefficients af ^/cos G appear to follow a single, power-law radius depen¬ 
dence (Fig. right panel). The coefficient aj j is the strongest outlier from this trend. For chip Nl, 
ajj « 1.0 X 10^^ whereas af q « 1.4 x 10^^, indicating that an additional physical mechanism acts 
on directly neighboring pixels.^ The moderate asymmetry of charge shifting is amplified by the 
differential nature of how the shift coefficients contribute to the pixel-to-pixel covariance to yield 
the strong asymmetry between Cov(l,0) and Cov(0,1) that we observed in Fig. ^ 


4. Effects on galaxy shape measurement 

One of the primary science applications of DES is to probe cosmic shear and the weak gravitational 
lensing effect of galaxies and clusters of galaxies. The measurement of gravitational shear y = 

^Compare also to the observed inverted resistivity dependence of the (0,1) lag. Note, however, that from the PTC 
flat series we find all covariances to be consistent with a g? scaling as predicted by the charge shift model. 
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12Ax + Ay 



Figure 8. Tests of the fitted shift coefficient model. Left panel; measured (black points with error bars) and 
model (red open squares) covariance coefficients for lags out to 5 for CCD N1. Right panel: comparison of 
slope of observed gain g w.r.t. flat count level /r to sum of model coefficients for all chips. According to 
Eqn. 3.8, we expect the two quantities to relate with a factor of 4g (solid red line). The mean agreement is 
better than at the ten per-cent level (dashed red lines). 





Figure 9. Shift coefficients of chip N1 plotted as function of Ax (left panel) and Ay (center). Projection 
effects are corrected in the right panel, where we normalize coefficients by the cosine of the angle between 
pixel border normal vector and direction of the lag. This leads to a well-behaved r^^-like power-law behavior 
(indicated by dotted green line) of the coefficients, with the neighbors and aj j (the latter is overplotted 
with the open symbol) as outliers. 


7i + iyi is based on a measurement of galaxy shapes e = ei+ i£ 2 - The latter can be defined based 
on second-order moments of the pre-seeing surface brightness profile 7(6) of a galaxy cenfered af 
0 = (O,O), 


jmeidjA^e 

.//(0)d2e ’ 


(4.1) 
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— Qn + Qii + '^x QnQii — Q 



(4.2) 


When second moments are measured with uniform weight as in Eqn. 0 deconvolution with the 
PSF can be performed by subtracting the observed and PSF second moments,^ 


Q 


,gaL dec _ 


-Q 


PSF 


(4.3) 


Observationally, one estimates 2*^^ from the second moments of star images. Since the stars 
used for this are typically at a high flux level in comparison to the galaxy images, the primary effect 
of charge-self interaction is a misrepresentation of the PSF second moments.^ 

In order to simulate the full effect on charge self-interaction on shape measurement, we imple¬ 
ment the A14 model in GalSim [30].^° Simulating PSF-convolved galaxy and star images with 


the effect applied, we estimate shapes using Eqns. |4. I| - p3| . We quantify biases in the measured 
shape relative to the input galaxy shape as 


gmeas_^true 

^meas_^true 


+ Cl 


+ C2 


-h 


p1 p? 
p\ pI 


(4.4) 


with multiplicative and additive biases m, and c, (e.g. |[l5[]) and leakage terms pi of the PSF 
ellipticity ef into the same and respective other component of shear. 

For the fiducial simulation settings, we choose the following conservatively realistic param¬ 
eters. The pre-seeing galaxy and PSF profile are modelled as Gaussians with full-width at half¬ 
maximum FWHMgai = 0.5" and FWHMpsp = 0.9", converted to the DECam pixel scale of 0.27" 
per pixel and Gaussian standard deviation a FWHM/2.35. We draw the star image with a 
FFUX_MAX in the central pixel of 60k charges. The galaxy is chosen to have a total flux of 4000 
charges. We add no noise to the image, assume a background level of Z? = 0 in the fiducial settings 
and set both galaxy and star image centroid to fall on a pixel center. For the shift coefficients, we 
use the mean values over the individually fitted chips out to lags A = 8pix. From determining 
at four different settings (where one of the components of either or is set to 0.1 and all 
others are 0) we fix the eight bias parameters of Eqn. 

Results for the fiducial settings and a number of variants are shown in Table ||. Multiplicative 
biases are above the per-cent level and exceed the DES weak lensing science requirements. For a 
hypothetical deep all-sky lensing survey, [|^ have calculated requirements on systematic errors in 


shape measurement that translate to approximate limits of \m\ > 0.1 x 10 


> 0.3 X 10 ^ and 


IpI > 0.3 X 10 all of which are exceeded by the observed biases. 

We probe the seeing range typical for DES lensing band data and object sizes between the 


smallest usable and median angular size of faint galaxies (cf., e.g., [128|]). We note that multiplicative 


“This is a simplified procedure for shape measurement, which in practice needs to use appropriate weighting to 
prevent the divergence of noise. Biases in real shape measurement pipelines should, however, be similar since the 
primary effect of charge self-interaction is a change in PSF size due to shifting of charges on scales of its FWHM. 

^Note that the effect can be alleviated by using fainter stars only, which are less affected by charge self-interaction 
(for DECam analyses, cf. e.g. [^). 

*®cf. https : //github . com/GalSim-developers/GalSim/ 
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bias is a strong function of the relative size of PSF and galaxy, as expected due to the PSF size 
misrepresentation in bright stars. High background levels somewhat alleviate the multiplicative 
bias since the shifting of background charges is a convolution effect that is independent of the 
object flux. Interestingly, the strong x/y asymmetry of neighbouring pixel covariances in flat fields 


(cf. Sec. 3.2) results in a relatively small additive bias c\ that disappears when we artificially set 


to the fitted value of 


-‘ 0.1 


1 . 0 - 


We apply reverse charge shifts to test the appropriateness of our correction scheme, e.g. use 
shift coefficients a —)• —a and apply Eqn. The purpose of this is to test the influence of a 
number of simplifications in the corrections we make on shape measurement biases (summarized 
in the second part of Table |I|). Specifically, 

• we correct according to the true model, but only out to lags of A^ax = 5pix; multiplicative 
bias is reduced by a factor of ps 30 to 8 x 10^^ for the fiducial case and is, as all other bias 
parameters, within the requirements; 

• we apply a maximally asymmetric model, degenerate in terms of flat field covariance mea¬ 
surements (aq j = 0 as transformed with Eqn. 3.10), but make a correction for the symmetric 


model fitted as described in Section |3.4| ; the most significant effect is that additive bias c\ is 
corrected only partially, but still below science requirements; 

we apply the model, using a Eanczos-3 interpolated version of the image at the pixel border 


centers, but for the reverse use the mean counts in the two neighbor pixels in Eqn. |3.2^ this 
leaves a per-mille level multiplicative bias, at the 5% level of the uncorrected case. 

We conclude that a A^ax = 5, symmetric, linearly interpolated correction scheme is sufficient for 
our scientific requirements. Eor a deep all-sky leasing survey, the latter two cause residual biases 
at the limit of acceptable systematic errors. 


5. Correction on pixel-level basis 

In the previous Sections ^ and ^ we have characterized flat-field noise pixel-to-pixel covariances 
with the A14 model and shown that the amplitude of charge self-interaction present in DEC am is 
expected to cause shape measurement biases at a critical level. Here, we test whether a correction 
using our model restores the profiles of bright stars to the shape of their fainter counterparts. 

To this end, we analyze a set of 90 r band exposures of moderately dense stellar fields taken 
between Aug 16, 2013 and Jan 19, 2014. To these we apply overscan and bias subtraction and a 
correction of inter-CCD cross-talk and signal chain non-linearity. Eor evaluation of the correction. 


we generate two versions of the frames, in one of which we apply Eqn. ^ with the negative of the 
fitted shift coefficients (e.g. a —)■ —a). Einally, we divide by a master flat field. 

On the reduced frames we run SEXTRACTOR 0 and split the catalog by brightness. We 
use the number of counts in the brightest object pixel, FLUX_MAX, with four bins centered on 
FLUX_MAX of 5000, 10000, 15000 and 20000 ADU, all sufficiently below the saturation level of 
DECam. We run PSEEx ||3|] on each of the catalogs to generate a model of the PSE as measured 
from stars of different flux. We remove a small number of frames for which the fitted PSE models 
have an irregular appearance visually. We finally stack the PSE models to get the mean PSE profile 


- 15 - 








settings 

0 

1 

Cl [10-3] 

p1[10-2] 


fiducial 

2.4 

-0.5 

-0.6 

-0.7 

FHWMpsf = 0.7" 

1.7 

-0.5 

-0.4 

-0.5 

FHWMpsf = 1.1" 

3.1 

-0.5 

-0.8 

-0.9 

FHWMgai = 0.3" 

6.9 

-1.4 

-1.7 

-2.0 

FHWMgai = 0.7" 

1.1 

-0.3 

-0.3 

-0.3 

2000 e^ background 

2.2 

-0.5 

-0.6 

-0.7 

symmetric ajj := a^i 

2.5 

0.0 

-0.6 

-0.7 

corrected out to A = 5 only 

0.0 

0.0 

0.0 

0.0 

corrected with flat-field degenerate model 

0.0 

-0.2 

0.0 

0.0 

corrected neglecting interpolation 

-0.1 

0.0 

-0.1 

0.0 


Table 1. Shape biases due to charge self-interaction in DECam as determined from image simulations. 
Multiplicative biases m, additive biases c and PSF leakage p, defined as in Eqn. 4.4, are measured with 
fiducial settings (see text in Section 4), and the described differences in parameters. Note that biases of 
approximately \m\ > 0.1 x 10 ^^, |c| > 0.3 x 10 ^^ and \p\ > 0.3 x 10 ^^ are problematic for a hypothetical 
deep all-sky lensing survey. Eor DBS, requirements are less stringent but clearly exceeded by the observed 
multiplicative bias. We always find m — nii Ri m 2 , C2 ~ 0 and pj ~ Ei ~ 0 omit these from the table. 


0.015 



Ax [pix] 



Ax [pix] 


0.000 


-0.075 ^ 


-0.090 


Figure 10. Difference image of normalized PSF as measured from stars at FLUX_MAX of 20000 and 5000 
ADU. Left; without correction, the brighter/fatter effect shows as a flux deficit in the central 9 pixels that 
is compensated by excess charge on an annulus at r: 2-4 pix radius. Right; after applying the charge shift 
correction. 


at each of the flux levels, ensuring that we only use frames where the number of stars in each 
FLUX_MAX bin is sufficient to fit a bilinear model of the spatial variation of the PSF that we 
always evaluate at the chip center position. 

Fig. shows a difference image of the FLUX_MAX « 5000 and 20000 star profile before and 


- 16 - 




















0.01 


2x10-8 

1.5x10-8 


CM 

<• 10-8 
w 

OxlO-’ 

0 

10“ 1.5x10“ 2x10“ 

peak flux [ADU] 

Figure 11. Relative change in flux in the central pixel of bright stars (left) and sum of squared residuals 
between bright and faint stars (right) before correction (red squares) and after applying the reverse charge 
shifts predicted from the model (black crosses). Both are plotted as a function of peak flux (FLUX_MAX) 
and measured relative to stars of FLUX_MAX Ri 5k ADU with (very small) bootstrapped error bars. 




0.004 


.X 0.002 
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-0.004 





Figure 12. Difference in radius (left) and ellipticity (center, right for the two components) between bright 
and faint stars before correction (red squares) and after applying the reverse charge shifts predicted from 
the model (black crosses). Quantities are based on second moments Qij measured with a Gaussian weight 
function of (7 = 3 pix width. Both are plotted as a function of peak flux (FLUX_MAX) and measured relative 
to stars of FLUX_MAX r; 5k ADU with bootstrapped error bars. 


after correction. The dominant brighter/fatter feature of a flux deficit in the central pixels is largely 
removed. 

We calculate a number of metrics to quantify how well the correction removes the flux depen¬ 
dence of the PSF. Metrics based on the reduction of the residuals in the individual pixels’ fluxes 
(see Fig. ITT] for flux deficit in central pixel and sum of squared deviations) show a reduction of the 
effect at the 90 per cent level. 

For lensing applications, however, the more relevant quantities are based on second moments 
of the PSF (cf. Eqns. 0ff). Figure |I^ shows residuals of three combinations of second mo¬ 
ments that correspond to the PSF size and the PSF ellipticity between brighter stars and stars of 
FLUX_MAX Ri 5000. The residuals at somewhat larger radii that remain in the corrected image 
(cf. Fig. 10) have a significant influence on second moments. While the ellipticity effects of 
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charge self-interaction appear corrected for, size residuals are overcorrected significantly. We have 
checked that this effect is not due to statistical or systematic differences in centroiding of the PS- 
FEx models. Shape biases associated with the observed residuals after correction are likely still 
problematic. 

We were not able to find the cause of the size residuals at the time of this writing. We hypoth¬ 
esize, however, that they could be connected to either a break-down of the model at low charge 
levels and large distances or one of the simplifying symmetry assumptions required in constraining 


the model parameters from flat field covariances (cf. Section 3.2). 


6. Summary 

We have presenfed a sysfemafic sfudy of imaging dafa from DECam and identified fwo effecfs of 
charge self-inferacfion: (1) a broadening of sfar images wifh increasing flux and (2) a correlation 
of Poisson noise in pairs of pixels in flal field exposures. 

The firsl effecl shows as missing flux in fhe cenfral pixels of brighf sfars fhaf instead appears 
af Ri3 pix disfance from fhe cenfer. The difference befween brighf and fainf sfars is approximately 
linear in flux and, af fixed fofal flux level, independenf of exposure time (Section 

We have measured fhe second effecl wifh high signal-lo-noise ratio ouf lo pixel-fo-pixel dis- 
fances of ?alO pix (Section ^). There is a « 5 : 1 asymmelry of correlations wifh the direct 
neighbor pixels on the same row or column. Correlation amplitudes for more distant pixels are 


almost isotropic with an approximate power-law fall-off 


OC f 


-2.5 


. Covariances differ from chip to 


chip at the « 20% level. The DECam chips with higher resistivity show larger correlations with 
their diagonal and distant neighbors and smaller correlations with the adjacent pixel perpendicu¬ 
lar to the read-out direction. There is a color dependence of the measured covariances, which is 
most severe for the latter pairs (more than 40% increase in correlation from g to Y band) but much 
weaker for other lags. 

Both effects are phenomenologically connected by the A14 charge shift model. Erom the latter 


measurements and a number of symmetry assumptions, we fit the model parameters (Section 3.3). 
We simulate the effect of the predicted charge shifts on shape measurement (Section Q) and find it 
to be significant for ongoing Dark Energy Survey weak tensing science, with multiplicative shape 
measurement biases above the per-cent level. 

Einally, we test the ability of the charge shift model to remove the flux dependence of the 
profiles of stellar images (Section We find that the flux dependence is strongly reduced in terms 
of the flux deficit in the central pixels and the overall sum of squared residuals. However, significant 
flux dependence of stellar sizes remains, indicating that our treatment of the charge self-interaction 
effects on inhomogeneous surface brightness images is yet incomplete. 
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